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Oh" General properties of anisotropic superconductors with mesoscopic phase separation are 

analysed. The main conclusions are as follows: Mesoscopic phase separation can be thermo- 



dynamically stable only in the presence of repulsive Coulomb interactions. Phase separation 
enables the appearance of superconductivity in a heterophase sample even if it were impossi- 
ble in pure-phase matter. Phase separation is crucial for the occurrence of superconductivity 
in bad conductors. Critical temperature for a mixture of pairing symmetries is higher than 



s 

q ■ the critical temperature related to any pure gap-wave symmetry of this mixture. In bad 
conductors, the critical temperature as a function of the superconductivity fraction has a 
bell shape. Phase separation makes the single-particle energy dispersion softer. For pla- 
nar structures phase separation suppresses d-wave superconductivity and enhances s-wave 
superconductivity. These features are in agreement with experiments for cuprates. 
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I. INTRODUCTION 



It is generally accepted that the majority of high-temperature superconductors, such 
as cuprates, possess two principal properties distinguishing them from the conventional 
low-temperature superconductors. These properties are mesoscopic phase separation and 
anisotropic pairing symmetry. 

Phase separation in superconductors implies that not the whole volume of a sample 
is actually superconducting but it is separated into regions of superconducting and normal 
phases. The latter may even be insulating. There exist numerous experiments confirming the 
occurrence of the phase separation in high-temperature superconductors, as is summarized 
in Refs. [1-4]. 

The phase separation is termed mesoscopic since the regions of coexisting phases form a 
kind of fog of clusters or droplets, whose typical sizes, corresponding to the coherence length 
l co h, are in between the mean interparticle distance a and the length of the sample L, so 
that 

a < l coh < L . 

These regions are intermixed, being randomly distributed in space. In general, the phase 
droplets are not static but can be dynamic, randomly fluctuating in time. In any case, 
whether they are static or not, their main features are the mesoscopic size and chaotic space 
location. Because of the random spatial distribution of the mesoscopic phase nuclei, they 
can be called heterophase fluctuations [5] . 

The mesoscopic phase separation is, actually, a very general phenomenon inherent to 
condensed matter [4-7]. This phenomenon happens in many systems, being responsible for 
a variety of unusual effects. For instance, it plays the key role in colossal magnetoresis- 
tant materials [8-10] and relaxor ferroelectrics [11,12]. Heterophase fluctuations, spatial or 
spatio-temporal, can exist in physical systems, without any noticeable external influence, 
thus, being self-organized [5,13]. The action of external forces can, of course, provoke the 
appearance of such mesoscopic fluctuations [5,14], making them more intensive. However, 
in general, the noticeable external perturbations are not compulsory, and heterophase fluc- 
tuations can really arise in a self-organized way. In some cases, these fluctuations can be 
triggered by infinitesimally small stochastic noise that always exists in all realistic systems, 
which are never completely isolated from their surrounding but are not more than quasi- 
isolated [5,15]. 

The possibility of mesoscopic phase separation in superconductors was advanced [16,17] 
yet before high-temperature superconductors were discovered [18]. Theoretical models con- 
firm that this phenomenon can be thermodynamically profitable, rending the heterophase 
material more stable [17,19-23]. 

Another specific feature of high-temperature cuprate superconductors is the anisotropy 
of the gap. A number of experiments point at the predominantly <i-wave symmetry of 
the superconducting order parameter [24,25], though in some cases one claims that the 
isotropic s-wave symmetry can be dominant [26-29]. The majority of experiments evidence 
the existence of the mixed s + d pairing in cuprates [25,30-38]. Several theoretical models, 
blending s-wave and ci-wave features, describe the s + d superconducting gap state and 
provide a reasonable explanation for various experiments [39-52]. Thus, the occurrence 
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of anisotropy in the gap of high-temperature cuprate superconductors seems to be well 
established. 

In the present paper, we suggest a model of superconductor, which combines two main 
features: mesoscopic phase separation and anisotropic pairing symmetry. We study inter- 
play between these characteristics. We do not narrow down the consideration by fitting 
parameters to a particular material, but we rather concentrate on the general properties of 
the model. The basic goal of the paper is to formulate the principal qualitative conclusions 
characterizing such a superconductor with both phase separation and an anisotropic gap. 

The treatment of the suggested model of superconductor is based on the theory of het- 
erophase materials [5] possessing the properties typical of the matter with mesoscopic phase 
separation. There are some important points that are worth emphasizing in order to better 
understand the following consideration. 

First of all, one should keep in mind that the basic spatial structure of matter is defined 
by ions forming a crystalline lattice. Charge carriers, such as electrons and holes, exist inside 
the given crystalline structure of a particular solid. Therefore the properties of all spatial 
characteristics, e.g. interaction potentials, are prescribed by a concrete crystalline structure 
of ions forming the lattice. 

Superconductivity or normal conductivity are the features of the charge carriers, reflect- 
ing the level of correlations between the latter. A solid with the same crystalline structure, 
as is well known, can be superconducting or not depending on the values of thermodynamic 
parameters. The occurrence of superconductivity of carriers does not substantially change 
the crystalline structure of ions. Thus, superconducting and normal phases of carriers may 
coexist inside the same crystalline lattice. 

The coexistence of different phases, typical of mesoscopic phase separation, means that 
the spatial regions of the sample are occupied by different phases, these regions, generally, 
have diverse shapes and random spatial locations. Then, in order to describe the properties 
of the sample as a whole, one has to average over phase configurations. The procedure of 
such a heterophase averaging is rather nontrivial, being analogous to the renormalizat ion- 
group technique, when one averages out one type of fluctuations with temporal or spatial 
scales that are distinct from another type of fluctuations. In our case, the heterophase fluc- 
tuations are mesoscopic, which distinguishes them from microscopic quantum fluctuations. 
All mathematical details of the heterophase averaging over configurations have been thor- 
oughly expounded in review [5]. In order that the reader could catch the main points of this 
procedure, these are sketched in the Appendix. 

After averaging over phase configurations, one obtains a renormalized Hamiltonian rep- 
resenting the phase replicas that would occupy the whole sample with a certain probability. 
In this way one comes to the picture where all characteristics do not involve anymore ran- 
dom spatial distributions but correspond to the averaged quantities resulting from their 
averaging over these random phase configurations. In the following sections, we deal with 
such averaged characteristics appearing after the heterophase averaging, whose essence is 
surveyed in the Appendix. 
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II. HETEROPHASE SUPERCONDUCTOR 



A superconductor with phase separation is a sample consisting of intermixed regions of 
different thermodynamic phases. Assume there are two phases, superconducting and normal, 
enumerated by the index v — 1,2. Let v = 1 correspond to the superconducting phase, while 
v — 2, to the normal phase. Quantum states of the phases pertain to the related spaces 7i u , 
which are the weighted Hilbert spaces [5]. The phases can be distinguished by their order 
parameters, as the gaps A u (k) in the momentum space, so that 

Ax(k)^0, A 2 (k) = 0. (1) 

Another way of distinguishing phases is by the associated order indices, which are defined 
for reduced density matrices [53] and have been generalized for arbitrary operators [54]. The 
order index for a bounded operator A is 



u{A) = 



log 


\\A\\ 
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Considering, in the place of A, p-particle density matrices p pu of the phases v — 1,2, we 
have the following [53,54]. For the superconducting phase, the order indices of odd density 
matrices are 

and those of even density matrices are 

u{p pl )= l - (p = 2,4,6,...). 

But for the normal phase, the order indices of all reduced density matrices are zero, 

^(P P2 ) = (p= 1,2,3,...). 

Coexisting phases occupy different spatial regions of the sample. These regions are 
composed of mesoscopic subregions that are randomly intermixed in space, forming com- 
plicated configurations. For each given configuration, we can define a locally-equilibrium 
Gibbs ensemble. Then, since the spatial phase distribution is random, it is necessary to 
average over these phase configurations. This procedure makes the basis of the theory of 
statistical systems with mesoscopic phase separation [5]. After averaging over random phase 
configurations, we come (see Appendix) to the renormalized Hamiltonian 

# = #x © # 2 (2) 

defined on the fiber space 

y = H 1 ®H 2 , (3) 

being the tensor product of the weighted Hilbert spaces. The phase-replica Hamiltonians 
H v can be written in the form 
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2 Tjint 



(4) 



where H* m is an operator of kinetic energy, if * n * is an operator describing pair interactions, 
and w u are phase probabilities satisfying the conditions 

Wi + w 2 — 1 , < w u < 1 . (5) 

The phase probabilities are defined as the minimizers of the thermodynamic potential 

Q — —T In Try e~^ k {(3T = 1) . (6) 

Here and in what follows, T is temperature, ks = 1- Setting the notation 

Wi = w , u>2 = 1 — w , (7) 

the minimization condition reads 

d 2 n 



<7U> 



> . 



(8) 



This condition shows when mesoscopic phase separation is profitable, as compared to a pure 
system. 

The average of an operator A is given as 



< A > = Try pA 



P 



Tr v e-f> H 



(9) 



with p being the statistical operator. Then the first equation from condition (8) takes the 
form 



dH 

< 7^- > = ' 
ow 



(10) 



and the second equation yields the inequality, being the condition of heterophase stability 

2 



< >-(]<: (— 
dw 2 \ dw 



> 



> 



(11) 



To illustrate the meaning of Eqs. (10) and (11), let us take the Hamiltonian (2) with the 
terms (4) and suppose that the kinetic part H km and the interaction part if * nt do not depend 
on w u . Let us define the mean kinetic energy K u and the mean interaction potential & u by 
the notation 



K u = < El in > , = 2 < Hl nt > 



int 



Then Eq. (10) gives the probability of the superconducting phase 

$ 2 + K 2 - K x 



w = 



$1 + $2 



(12) 



(13) 



From here, since < w < 1, we get 
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-$i < K x - K 2 < $ 2 



(14) 



The stability condition (11) yields 

from where it follows that the necessary condition for the stability of a phase-separated 
sample is 

$i + $2>0. (15) 

In that case, phase separation becomes thermodynamically profitable, for which, as is seen 
from Eq. (15), the existence of repulsive interactions is compulsory. 



III. STRUCTURE OF HAMILTON! AN 

Employing the field representation, we deal with the field operators tp sl/ (r), in which 
s =j, | denotes spin and r is a spatial vector. Fermi commutation relations are assumed. 
The kinetic part has the standard form 

H kin = £ I ^ (r) ^ (r) _ ^ ^ (r) dr ) (16) 
s J 

where K v {y) is a kinetic transport operator and //, chemical potential. The interaction part 
H™*, in general, consists of direct interactions and of effective interactions due to a kind 
of boson exchange. This, for instance, can be the phonon exchange, if one considers the 
usual picture based on the Frolich Hamiltonian [55-58]. In principle, one may consider the 
exchange by other types of bosons, say excitons, but for concreteness, we shall keep in mind 
the conventional phonon picture. For simplicity, and at the same time for generality, we 
take for the interaction Hamiltonian the expression 

H l nt = \ E / C(r)Vi,(r')K(r, r')W„(r')?Mr) dvdr' , (17) 

ss' 

where the vertex operator V v models all effective interactions, direct as well as those caused 
by boson exchange. The vertex operator is supposed to be symmetric, 

K(r,r') =K(r',r) . (18) 

In agreement with condition (1), the anomalous averages for the superconducting phase 
are not trivial 

< V s i(r)Vvi(r') > ^ , (19) 
at least for some spins, while such averages for the normal phase are identically zero, 

< VWrM^r') > = 0. (20) 
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For crystalline matter with a periodic structure, the field operator can be expanded over 
Bloch functions, which, for the single-zone case, writes 

^su(r) =^2c av (k)(p k (r) , (21) 

k 

with k being wave vector. Let us introduce the matrix elements over the Bloch functions 
V?fc(r), resulting in the transport matrix 

t v (k,p) = (<p k ,k v <p p ) (22) 

and the vertex 

K(k, k', p', p) = ((p k (pk', KvVV^p) . (23) 
The latter, due to Eq. (18), has the symmetry property 

K(k,k',p / ,p) = K(k',k,p,p / ). (24) 
Invoking expansion (21), the kinetic part (16) transforms to 

H kin = £ £ UK p) _ ^ 4( k ) Cw ( p ) (25) 
s fcp 

and the interaction term(17) becomes 

C'^EEE ^( k > k', p', p)cL(k)cL(k')c^( P ')c s ,( P ) . (26) 

ss' fefc' pp' 

To make the problem treatable, let us resort to the Hartree-Fock-Bogolubov approxima- 
tion, according to which the four-operator products are expressed as 

ClC 2 C 3 C 4 = CiC 2 < C3C4 > + < CiC 2 > C3C4- < CiC 2 X C 3 C4 > + 

+C1C4 < c 2 c 3 > + < C1C4 > c 2 c 3 - 

- < C1C4 >< c 2 c 3 > -C1C3 < c 2 c 4 > - < C1C3 > c 2 c 4 + < C1C3 >< c 2 c 4 > , (27) 

where q represents any of the operators c su (k) or cj^(k) and the Fermi commutation relations 
are assumed. Also, we shall consider the restricted spaces of quantum states, for which the 
Bardeen-Cooper-Schrieffer restriction is valid: 

c L( k ) c sv(k') = S sg/ <Wc4(k)c S!y (k) , 

4(k)cL(k') = 5^ 5_ w 4(k)ci w (-k) . (28) 

This means that the restricted spaces consist of the wave functions for which spin and 
momentum are conserved. 
The normal average 
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ra„(k) = ]T < c L(k)c s ,(k) > 



(29) 



is the momentum distribution of particles. Introducing the anomalous average 

<7„(k) = < c_ w (-k)c w (k) > , (30) 
conditions (19) and (20) can be rewritten as 

<7!(k)^0, <7 2 (k) = 0. (31) 

With the approximations (27) and (28), the Hamiltonian (4) can be diagonalized by 
means of the Bogolubov canonical transformation 



c s „(k) = u v (k)a ai/ (k) + u v (k)aL ai/ (k) , 



in which 



|u„(k) 



1 + 



^(k) 
E u (k) 



K(k) 



1 - 



^(k) 
E v ik) 



Here, the single-particle dispersion is 

w„(k) = t„(k, k) + w v M v {k) - n , 

with the mass operator 



M v (k) = E K(k,p, P ,k)- -K(k,p,k,p) 
P L 1 



Mp) > 



and the excitation spectrum 



contains the gap 



A„(k) =w L ,J2 J »(k, p) ov(p) , 



(32) 



(33) 



(34) 



(35) 



(36) 



where the effective interaction 

J u (k, p) = -K(k, -k, -p, p) . 
Then the Hamiltonian (4) reduces to 

H v = w u ^ ^ZE v {k)a\ u (k) a su{k) + w v C v , 

s k 



with the nonoperator term 



uj v {k) - E u ik) + A,(kK(k) - - w u M u {k)Mk) 



(37) 



(38) 



For the averages (29) and (30), one gets 



/, \ i ^( k ) , i w v E„(k) A„(k) ^^(k) 

n,(k) = 1 " ^ tanh , a,(k) = _ tanh . (39) 

These expressions have sense for both the superconducting and normal phases; however, for 
the normal phase, according to condition (31), one has 

<r 2 (k) = A 2 (k) = , (40) 

thence 

2 

n2{k) = exp{(3w 2 u 2 (k)} + 1 ' 

Note that the phase probabilities w u enter all equations in a rather nontrivial way, which 
will essentially influence the properties of superconductors with mesoscopic phase separation. 



IV. PHASE SEPARATION 

In what follows, we shall use the notation (7), writing w — w±, and, similarly, we shall 
omit, for the sake of simplicity, the index v — 1 at all related quantities. For instance, we 
shall write A(k), <x(k), E(k), and so on instead of Ai(k), <7i(k), and -E'i(k). 

The gap equation (36) can be presented as 

A(k) = ^£j(k, P )^tanh^M. (41) 

Looking for a positive solution for A(k), we see that this is possible when the right-hand 
side of Eq. (41) is also positive, which requires that the effective interaction 

J(k,p)>0 (42) 

be positive in the region of momenta making the main contribution in the summation of Eq. 
(41). 

Another necessary condition is the condition (15) for the profitability of phase separation. 
For the mean interaction potential, defined in Eq. (12), we find 

i = ^M,(kK(k)-2^J„(k )P K(kK(p) . (43) 

k kp 

Then the stability condition (15) yields 

]T £ M w (k)n„(k) > 2 ]T J(k, pMk)a(p) . (44) 

v k kp 

From here, we get the necessary condition 

£ ]Tii4(kK(k)>o, (45) 

v k 

which tells that some sufficiently strong repulsive interactions are to be present in the system. 
Such natural interactions are, of course, Coulomb interactions. Thus, we come to the first 
conclusion: 
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Mesoscopic phase separation in superconductors can be thermodynamically stable only in 
the presence of repulsive Coulomb interactions. 

Now let us recall that, as is discussed in the Introduction, the gap of the hole-doped 
high-temperature cuprate superconductors displays strong anisotropic dependence on mo- 
mentum. To describe the anisotropy, one may introduce a basis {x«(k)} of functions Xi(k) 
characterizing the lattice symmetry, with the index i — 1, 2, . . . enumerating irreducible rep- 
resentations of the symmetry group. Let such a basis be defined, being orthonormal and 
complete, 

E X*(k)Xj(k) = , Yl X*(k)Xi(p) = hp ■ 



Then one may expand (see [25], [59]) over this basis the effective interaction 

J(k,p) = ]Tj^(k) x *(p) 



V) 



and the gap 



A(k) = EA*(k). 



(46) 



(47) 



Using this, the gap equation (41) reduces to 



(48) 



where 



= W r f, t \ tanh 
3 p 2£(p) 



wE(p) 
2T 



x*(p)xj(p) • 



(49) 



The system of uniform algebraic equations (48) possesses nontrivial solutions when 

det(l-i)=0, (50) 

where 1 = [8^] is the unity matrix and the matrix A = [A^] is composed of elements (49). 

The effective interaction (46) consists of an attractive part, caused by phonon exchange, 
and a repulsive part, due to direct Coulomb interactions [60], because of which Jij has the 
structure 



J i/ 



\a\ 



M h 



(51) 



in which a is the charge-lattice coupling, c^o is the characteristic lattice frequency in the 
presence of heterostructural fluctuations [61,62], connected by the relation 



(52) 



with the characteristic lattice frequency c^o of a pure sample, and M is an effective intensity 
of direct Coulomb interactions. The latter approximately equals 



M 



7re , 
•,2 111 



1+4 t£ 



K 
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where ^ is a Fermi momentum of charge carriers, k" 1 is a screening radius, for which 
k 2 m 4moe 2 (3p/"7r) 1 / 3 , and m , e, and p are mass, charge, and density of carriers. 

Keeping in mind inequality (42), we set Jjj > and 6^ > 0. Then the condition for the 
existence of superconductivity reads 



\a 



, ,2 
uj 



2 

- wM > . (53) 



This differs from the Bardeen-Cooper-Schriefer criterion for superconductivity [63] by the 
presence of the superconducting phase probability w, which essentially changes the meaning 
of Eq. (53). It may happen that Coulomb interactions are so strong, with M > \a/u \ 2 , 
that superconductivity in a pure sample is impossible. However, since w < 1, criterion (53) 
may be valid, which implies the occurrence of superconductivity. In this way, we get the 
second conclusion: 

Phase separation enables the appearance of superconductivity in a heterophase sample 
even if it were impossible in pure-phase matter. 

Defining the dimensionless quantity 

V* = M ^| , (54) 



\a\ 



condition (53) can be presented as 

1 - wp* > . (55) 

The parameter (54) is of the order [60] of p* ~ ujq/uj p , where u p is the ion plasma frequency. 
For good conductors, c^o *C iv p , hence /i* C 1, and inequality (55) is easy to satisfy even for 
a pure sample, with w — 1. For bad conductors, c^o > uj p , so that pf > 1. In particular, 
if pf > 1, superconductivity cannot arise in a pure sample, though may appear in phase- 
separated matter, with w < 1. This yields the third conclusion: 

Phase separation is crucial for the occurrence of superconductivity in bad conductors. 

V. CRITICAL TEMPERATURE 

At the critical temperature T c , the gap tends to zero, A(k) — > 0, hence E(k) — > u;(k). 
Then Eq. (49) becomes 



AijfTc) = V^^tanh 



WLd(p) 

2T C 



Using the density of states 

JV«M = E% - w(p))x,-(p)%(P) . ( 57 ) 
p 

satisfying the normalization 

r+oo 

Nij(u) du = 5ij , 



11 



Eq. (56) can be written as 

MT C ) = wJ l3 £J M± tanh du . (58) 

The density of states (57), with the standard replacement of summation by integration 

v - f dp 

where B implies the Brillouin zone, transforms to 

= -J B S(u - W (p)) XI(P) XAP) ^ ■ (59) 

Assuming, as usual, that the density of states N^u) is the largest on the Fermi surface and 
fastly decreases after uj > u , Eq. (58) can be reduced to 

Aij(T c ) = wJ tJ N tJ (0) / - tanh — ) du . (60) 
Jo uj \Z1 C J 

Introducing the coupling matrix A, with the elements 

l«l 2 

Xij = N^oy-L hj , (61) 

effective coupling matrix A, with 

A,,- = wJijNijiO) = (1 - wn*)^ , (62) 
and the characteristic integral 

/c " X x tanh i^r x ) dx ' (63) 

we can present Eq. (60) as 

A l3 (T c ) = I C A^ ■ (64) 
Substituting this into condition (50) gives the equation 

det(l - I C A) = (65) 

for the critical temperature T c . 

If Ay is diagonal, then T c is defined by the largest Ajj. However, in general, for anisotropic 
superconductors, A^ is not diagonal, and A^ ^ for % ^ j. The latter means that, generally, 
the gap (47) is presented by a mixture of waves of different symmetry. There is a very 
nontrivial relation between the gap being such a mixture and the magnitude of the critical 
temperature, which we describe below. 

Let us define an effective coupling A e ff by the identity 



12 



A eff I c = 1 - det(l - I C A) 



(66) 



Then Eq. (65) for the critical temperature takes the form 

Kfflc = 1 • (67) 

Note that, since, according to definition (63), the integral I c > is positive, then A e ff > 0. 
From Eq. (67) it follows 

8T IT 2 1 

' " 77 > " • (68) 



6A e ff LUoW Z l 2 A eff 

where 



J' = / sech 2 ( - 1 2 r^ Q x) dx . 



Inequality (68) tells us that T c is higher for larger A e ff. This is valid for all A e ff > and 
can be explicitly illustrated for the particular cases: 

T c ~ 1.14 w 3 ' 2 u Q exp (- ^—^j (A eff < 1) , 

J w 3 / 2 uj A eff (A eff » 1) . (69) 



" 2 

Another important inequality is 

A eff > maxAji , (70) 

i 

provided that A^- 7^ for some i 7^ j. And, if A^- = SijAu, then A e ff = maxjAjj. This 
property is easy to explicitly demonstrate for the case, when there are two prevailing wave 
symmetries. Thus, if i = 1,2, then Eq. (66) yields 



1 

Ae// = g 



A 11 + A 22 + a/(Ah-A 2 2) 2 + 4A 



12 



(71) 



Without the loss of generality, we may use the enumeration such that An > A 22 - As is seen 
fromEq. (71), 

Aeff = A n (A 12 = 0) . 

In addition, 



dA 



eff 



dA\ 2 ,/(A 11 -A 22 )2 + 4A? 2 



> , (72) 



hence A e ff increases with increasing |A 12 |. Thus, for a mixture of waves, when |A i3 -| > 0, 
where % 7^ j, the effective coupling A e ff becomes larger than the maximal A u . But, in 
agreement with Eq. (68), the larger A e jj, the higher is the transition temperature T c . 
Therefore, we come to the following conclusion: 



13 



Critical temperature for a mixture of gap waves is higher than the critical temperature 
related to any pure gap wave from this mixture. 

In superconductors with phase separation, the critical temperature can be a nonmono- 
tonic function of the superconducting fraction w. To show this, let us consider the case 
when one of the gap symmetries is prevailing, so that one of is essentially larger than 
other Ajj. Let us denote this maximal An as 

maxAjj = (1 — wfi*)X , (73) 

i 

where the relation (62) is taken into account. The equation (67) for the critical temperature 
can be written as 

(1 - Wfl*)\I c = 1 . 

The characteristic integral (63) possesses the following asymptotic properties: 

I.*^ <«^0), 

/^l„(l.l4^) (JWO). 

Using these properties, we see that the critical temperature tends to zero in two limits, if 
the superconducting fraction tends to zero, when 

T r 1 



■ C 



[1 - wtfXw 3 * 2 («;->0), (74) 
and also if this fraction tends to a finite value 1///*, when 

— ~ 1.14iw 3/2 exp|- — r-1 [iu->— ]. (75) 

For good conductors, when /j*<1, the limit (75) is unachievable. But for bad conductors, 
for which /i* > 1, this limit can be achieved. This can be formulated as another conclusion: 

In bad conductors, the critical temperature as a function of the superconducting fraction 
w has the bell shape, tending to zero at w — > and at w — > . 

To estimate the point where the critical temperature is maximal, we may keep in mind 
that experiments with high-temperature superconductors show that only part of a given 
sample is in a superconducting phase, this part often being just a few percent [64,65]. This 
means that w < 1 and, hence, the value w max , where T c = T max is maximal, is also small, 
Wmax 1. Taking this into account, from the above equations we obtain 

Tmax ^ - A UJ W^ X , W max ~ — . 

Such a bell shape of the critical temperature as a function of doping is typical of experi- 
mental curves for high-temperature cuprate superconductors [1, 66,67], where the maximal 
critical temperature occurs at the optimal doping 0.15. Assuming that the superconducting 
fraction is proportional to the doping, so that w max ~ 0.15, we have /i* ~ 4. Then the func- 
tion T c (w) has a striking similarity with the behaviour of T c as a function of doping, studied 
in experiments with cuprates. Figure 1 illustrates T c (w) found by solving numerically Eq. 
(67). 
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VI. DENSITY OF STATES 



To concretize the consideration, let us take into account that the crystalline structure of 
cuprates is such that the carriers move mainly in planes, only rarely jumping between the 
latter, which are separated by a distance essentially exceeding the mean distance a between 
lattice sites on the plane. Neglecting the interplane jumps, one comes to a two-dimensional 
motion of carriers on the plane. In the case of such a planar motion, the single-particle 
dispersion cu(p) = ^i(p), given by Eq. (33), depends only on two momentum components, 
say pi and p 2 . Then in the following formulas, it is easy to integrate out the third component 

In many cuprates, the Cu and O atoms arrange themselves in a square lattice with the 
point group symmetry C^ v [25]. This concrete case will be employed in what follows. 

For a square lattice a x a, the Brillouin zone is defined by the wave vectors p a G 
[— n /a, ir/a], with a — 1,2. It is convenient to introduce the dimensionless wave vectors 
k = {k a }, in which k a = p a a, so that k a G [— 7r, 7r]. The point group C 4l) of a square lattice 
is characterized by three one-dimensional irreducible representations labelled as Ai, Bi, 
and B 2 . The representation A\ is of rank 3, having three types of symmetries denoted as 
s, s*, and s xy . The representation B\ is of rank 1, with the symmetry d x 2_ y 2. And the 
irreducible representation B 2 is of rank 1, with the type of symmetry denoted by d xy . The 
corresponding basis functions are 



Xi(k) = 1 , (s) 

X 2 (k) = ^ (cosfci + cosk 2 ) , (s*) 

X 3 (k) = cos h • cos k 2 , (s xy ) 

X 4 (k) = sin fci • sin/c 2 , (d xy ) 

Xs(k) = ]- (cosfci - cos^) . (d x 2_ y 2) 



z 

For the density of states (57), we have 

JVyM = £ 8(f* ~ ^(k))x*(k)x,(k) . (76) 

The basis functions xl(k) = X«(k) are real and symmetric with respect to the inversion of 
k, so that Xi(~ k) = Xi(k). The single-particle dispersion is also symmetric, u(— k) = u;(k). 
It is convenient to introduce the dimensionless dispersion 

S (k) . ^ . (77) 

Then the density of states (76) on the Fermi surface becomes 

%(°) = ^ £ *(w(k))x<(k)Xi(k) dfci ^2 • (78) 
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To write down an explicit expression for the single-particle dispersion w(k), one usually 
resorts to the tight-binding approximation [25,39,40,45], in which for a square lattice, one 
has 

uj(k) = — t e ff (cos ki + cos k 2 ) — [i , (79) 

where t e ff is an effective transport parameter. In our case, taking into account Eq. (33), 
we see that t e jf consists of two parts, 

t eff = t + wM , (80) 

the nearest-neighbour hopping integral t and the intensity of the repulsive Coulomb inter- 
action M . These terms correspond to the transport matrix and mass operator, respectively. 
As is well established, both experimentally and theoretically, a strong Coulomb repulsion is 
present in all cuprates [25]. With the parameter 

t = hfL = i(t + W M ), (81) 

the dimensionless dispersion (77) takes the form 

57(k) = -*(cos h + cos k 2 ) - 1 . (82) 

When the mesoscopic phase separation happens in a superconductor, so that w < 1, then, 
as is seen from Eq. (80), the parameter t e ff decreases. As a result of this, the dispersion 
(79) becomes softer. Thus, we get an important conclusion: 

Phase separation softens the single-particle dispersion. 

To proceed further, we calculate the density of states (78). For this purpose, we use the 
dispersion (82) and note that 

v v ;; |i sin A; 2 (A;i)| 
where the function k 2 (k) is defined by the equation 

cos k 2 (k) = — - — cos k . 

t 

Then the density of states (78) can be transformed to 

nm = 4- r ^ k ^ dk , (83 ) 

^ Jk _ (! + t cc .s A;) 2 

where ^ 

fi(k) = Xi(k, k 2 (k)) , k = arccos ( 1 - - 



< /,•() < 7T ( i < t < OO 
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Explicit expressions for the functions (fii(k) are 

1 



(pi(k) = 1 , <p 2 (k) = 



2t ' 



+ cos k ) cos k 



ifi(k) = sin k y 1 — + cos fc^ , 
Accomplishing in Eq. (83) the change of variables 

ipi(x) = ipi(&rccosx) (x = cos/c) , 

we obtain the form 



¥h{k) = — + cos A; 



%(0) 



TT 2 jJ,t 



dx , 



(84) 



in which 



^i(x) = 1 , ^(x) = - ^ > ^a(^) = - Q + ar) a; 
V>4(z) = 



\ 







(1 - a: 2 ) 


1 -G +I )' 







The value of the quantity iVy(O) plays an important role in defining the coupling parameters 
(61) and (62). The larger is the density of states Na(0), the more profitable is the occurrence 
of the related gap symmetry labelled by the index i. 

Let us analyse the behaviour of Eq. (84) as a function of t changing from t = 1/2 to 
larger t > 1/2. For convenience, we consider the dimensionless quantity 



Dij = 7r fj,Nij(0) , 



(85) 



normalized by means of 



N n (0) = N 22 (0) 



TTfl 



t 



1 



At t — 1/2, integral (84) can be calculated analytically, giving 

D u = D 13 = D 22 = D 33 = 1, D 14 = D 24 = D 34 = D M = D 55 = 

V 



Dis — D 22 — D 33 — 1 , 
D\2 = D 23 



Also, for alH > 1/2, one has 



D 



15 



D 



25 



-1 



D 



t 



35 



D. 



45 



. 



For t > 1/2, we calculated the integral (84) numerically. The corresponding densities of 
states are shown in Fig. 2. Among all D u , only D n and D 22 monotonically decrease, all 
other Da are nonmonotonic functions of t. Comparing the density of states Da for different 
% = 1,2,3,4,5, we have the following. Since at t — 1/2, only D n = D 22 = D 33 = 1 are 
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nonzero, while D44 = D55 = 0, thence only the s, s*, and s xy waves can exist. But at t — 1, 
we get D n ks 0.6, D 22 = D44 ~ 0.2, D 33 ps 0, _D 55 ps 0.1, which tells that now the gap 
symmetries s, s*, d xy , and d x 2_ y 2 may exist, while the s xy wave disappears. With increasing 
t, the density D 55 increases. For instance, at t — 3, we have -D n pa 0.3, D 22 , and _D 33 are 
close to zero, D44 = D 55 Pi 0.1. Therefore, here the probable symmetries are s, d xy , and 
d x 2_ y 2. For t > 3, the highest densities are Dn, diminishing below 0.3, and -D55 ~ 0.1, all 
other densities being smaller. In this way, for t > 3 the most probable symmetries are s and 

d x 2— y 2 . 

To estimate the magnitude of t for high-temperature superconductors, we may take the 
values of parameters typical of cuprates [25,39,40,45], that is, t ps (0.5 — 1) eV, M ps 1 eV, 
and n ~ 0.5 eV. Then the parameter t, defined in Eq. (81), with w ~ 1, equals t ~ 3 — 4. In 
this region of t, the largest densities are Dn and -D55, hence the symmetries s and d x 2_ y 2 are 
preferable. But when phase separation occurs, the superconducting fraction w becomes less 
than unity, which diminishes the value of t. The decrease of t suppresses the density of states 
D 55 and enhances Dn, which means that the relative weights of s and d x 2_ y 2 symmetries 
are changing. This can be formulated as the following conclusion: 

Mesoscopic phase separation suppresses the contribution of d-wave superconductivity and 
enhances that of s -wave superconductivity. 

The value of the superconducting fraction w depends on such parameters as temperature, 
pressure, or external magnetic fields. Therefore the relative contribution of different wave 
symmetries will be varying under the action of these parameters. This can explain why in 
different experiments one observes alternately the dominance of either s or d gap symmetries. 

VII. SUMMARY 

We have studied the main properties of superconductors, such as cuprates, exhibiting two 
principal features common for these high-temperature superconductors, mesoscopic phase 
separation and anisotropic gap symmetry. Interplay between these two phenomena is in- 
vestigated by means of a model incorporating the presence of mesoscopic phase separation 
into the randomly distributed regions of superconducting and normal phases. The following 
general conclusions are obtained: 

1. Mesoscopic phase separation in superconductors can be thermodynamically stable only 
in the presence of repulsive Coulomb interactions. 

2. Phase separation enables the appearance of superconductivity in a heterophase sample 
even if it were impossible in pure-phase matter. 

3. Phase separation is crucial for the occurrence of superconductivity in bad conductors. 

4. Critical temperature for a mixture of gap waves is higher than the critical temperature 
related to any pure gap wave from this mixture. 

5. In bad conductors, the critical temperature as a function of the superconducting frac- 
tion has the bell shape. 
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6. Phase separation softens the single-particle energy dispersion. 

7. Mesoscopic phase separation suppresses the contribution of d-wave superconductivity 
and enhances that of s-wave superconductivity. 

These conclusions are in good qualitative agreement with experiments for high- 
temperature superconductors. Since in colossal magnetoresistance materials there also oc- 
curs the phenomenon of mesoscopic phase separation [8-10,68,69], such materials may pos- 
sess some of the features described in this paper. Moreover, the mesoscopic phase separation 
is a rather general phenomenon appearing in different kinds of condensed matter and it can 
be described in the frame of the general theory [5,70], which was used in this paper for 
studying the general properties of anisotropic phase-separated superconductors. 
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Appendix. Averaging over Phase Configurations 



Let the considered sample occupy in the real space {r} a region V having the volume 
V = fydr. Assume that mesoscopic phase separation occurs in the sample, so that it is 
filled by several thermodynamic phases, which are enumerated by an index v. Say, v — 1 
corresponds to superconducting phase, while v — 2, to normal phase. At each given instant 
of time, distinct phases are located at different spatial regions, which can also be labelled by 
the phase index v. This means that the total sample volume V can be divided into subregions 
V„, filled by the related phases. A family {Y u } of subregions Vj, forms an orthogonal covering 
of V. This covering can be characterized by a family 

£ = {ar)|reV} 

of the manifold indicator functions 



1, reY u 
0, r £ V, . 



The family £ uniquely defines a phase configuration in the real-space volume V. 

From the physical point of view, distinct thermodynamic phases possess different prop- 
erties, because of which such phases can be distinguished from each other. For example, 
the phases can be distinguished by their order parameters. In the case of superconductor, a 
convenient order parameter can be chosen as the anomalous average 

r) u (r) = < ^(r)^(r) > £„(r) , 

where ip v (r) is a field operator of the field of carriers, when the latter are in the phase 
u, and £„(r) is the related indicator function. For the superconducting regions, one has 
r)i(r) 7^ 0, provided that r 6 Vi, while for the normal parts, one has 772(f) = 0, with 
r G V 2 . The anomalous average, as is known, is directly related to the gap in the spectrum 
of excitations. Thus, different phases could be distinguished by the existence or absence of 
the gap in the spectrum. For a while, when dealing with a nonuniform sample consisting 
of separate phases, the anomalous average, depending on the spatial variable r, is a more 
convenient order parameter. 

Mesoscopic phase separation occurs in the real space in a random way, which means that 
the locations and shapes of the phase subregions Y u are random. Then observable quantities 
should be determined with averaging over these random phase configurations. Since each 
phase configuration is uniquely defined by the set £ of the manifold indicator functions, it 
is necessary to describe an ensemble {£} of all possible sets £, corresponding to all possible 
phase configurations. For this purpose, we introduce an orthogonal subcovering {V^} of 
each region V„, such that 

™ v — u i=l v vi 1 v fit 1 1 * vj — <Jfiv u t] v vi ■ 

To each subregion V^, we ascribe a vector a ui e V„i, called the center, playing the role of 
a local center of coordinates, so that moving a ui implies a congruent motion of V„j. The 
introduced subcovering is uniquely characterized by a family of the indicator functions 
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with the property 



1, reV* 
0, r £ V ui , 



i=i 

By moving the centers and changing the measure of V„, it is possible to construct various 
phase configurations. To explicitly realize the averaging over these configurations, we define 
the differential functional measure 

in which 

x v = — / f„(r) dr . 
V Jv 

This measure, with varying a ui G V and x u G [0,1], induces a topology on the manifold {£}, 
which results in the topological configuration space X = {£| composed of all admissible 
phase configurations. 

For each fixed phase configuration, the observable quantities, represented by Hermitian 
operators, depend on the given configuration and have the structure 

ah) = ®u A v {i u ) , 

being defined on the fiber space y = ® V H. V , whose fibering yields the weighted Hilbert spaces 
H. v . The dependence of the operators of observable quantities on the indicator functions, 
marking the space filled by the corresponding phases, naturally comes from the identity 



/ dr = [ £„(r) dr . 



This identity is employed when representing the operators through the integrals of their 
operator densities: 

i(0 = / A(£,t) dr = ® v A„(S V ) 

A(Z, r) = @ v Au(tu, r) , A v (^ v ) = J A u (£ u , r) dr . 

Here and everywhere in the paper, we denote, for simplicity, the integration over the whole 
sample as 

For example, let us write down the Hamiltonian density 

H v fo,r) = UrM(r)Ku(r)Mr) + \ J ^(r)6(r , )^(r)^(r')K(r, r')^(r')^(r) dr dr' , 

in which K v {r) is an operator of kinetic energy and V„(t, r') is an effective interaction. 
Analogously, the number-of-particle operator density is 

iV,(6,r)=6(r)^(r)^(r). 
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The field operators here are assumed to be columns with respect to spin indices. 

A nonuniform system, composed of several thermodynamic phases, must be described 
by the quasiequilibrium (or locally equilibrium) Gibbs ensemble, with a statistical operator 
proportional to e~ x ^\ where 



X(0 = J r) [H(t, r) - r)JV(£ r) 



dr 



Here the local inverse temperature r) and the local chemical potential r) model the 
system nonuniformity corresponding to a given phase configuration characterized by a set 
£. The natural thermodynamic potential for a quasiequilibrium system, with random £, is 



Q = - In Ti y J e~*® . 



Assuming, as usual, the existence of the thermodynamic limit, and accomplishing the aver- 
aging over phase configurations, which is characterized by the differential measure D£, it is 
possible to prove [5] that the thermodynamic potential Q reduces to the form 

Q = - In Ti y e~ p6 , 
with an effective renormalized Hamiltonian H = (B U H U , in which 

H v = / H y {w u , r) - {jlN v {w v , r) dr , 

and where the average inverse temperature and average chemical potential, respectively, are 

(3 = fp(t,r)Dt, ii = JiM(t,r)Dt. 

Then the potential 

Q = (3tt 

is simply related to the Gibbs grand potential Q. The phase probabilities w u are defined as 
the minimizers of either Q or Q. 

In this way, after averaging over heterophase configurations, we come to the description 
of the system by means of a renormalized Hamiltonian, containing only averaged quanti- 
ties and not involving anymore random phase distributions that have been averaged out. 
All expressions throughout the paper correspond to averaged quantities resulting from the 
described procedure of heterophase averaging. A complete and detailed mathematical foun- 
dation for this averaging procedure is given in reviews [5,70]. 
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Figure Captions 



Fig. 1. Critical temperature T c , in units of cj , as a function of w for //* = 4 and 
different couplings: A = 1 (dotted line), A = 5 (dashed line), and A = 10 (solid line). 



Fig. 2. Density of states Da as a function of the effective transport parameter t: Du 
(upper dashed-double-dotted line), D 2 2 (dashed-dotted line), D 33 (dotted line), D i4 (dashed 
line), and D 55 (solid line). 
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